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Alfred Stadler^'^* and Franz Gross^'^ 
^ College of William and Mary, Williamsburg, VA 23185 
^Centra de Fisica Nuclear da Universidade de Lisboa, 1699 Lisboa Codex, Portugal 
^ Thomas Jefferson National Accelerator Facility, Newport News, VA 23606 

(February 9, 2008) 

First results for the triton binding energy obtained from the relativistic spectator or Gross 
equation are reported. The Dirac structure of the nucleons is taken into account. Numerical results 
are presented for a family of realistic OBE models with off-shell scalar couplings. It is shown that 
these off-shell couplings improve both the fits to the two-body data and the predictions for the 
binding energy. 

21.10.Dr,21.45.+v,21.30.-x,13.75.C,11.10.St 



Preprint # WM-96-109, TJNAF-TH-96-13, CFNUL-96-5 



1 



The first realistic nonrelativistic calculations of the 
triton binding energy were completed in the 1970's [Q. 
Later it was shown that different methods gave the same 
results, and that the binding energy could be calculated 
to an accuracy of a few keV by considering all nucleon- 
nucleon (TV TV) partial waves up to j = 4 [Q. Today, 
if three-body forces (3BFs) are not considered, a small 
discrepancy of about 0.5-1.0 MeV remains between the 
experimentally observed value of —8.48 MeV and val- 
ues obtained from realistic nonrelativistic NN potentials. 
State-of-the art calculations now include sophisticated 
3BFs, and when their strength is adjusted to give the 
correct triton binding energy, an excellent value is also 
obtained for the "^He binding energy (and to a lesser ex- 
tent other light nuclei up to A ~ 7) [||. 

However, relativistic effects should make a contribu- 
tion to the binding energy at the level of several hundred 
kcV. Using a mean momentum of about 200 MeV (con- 
sistent with nonrelativistic estimates) we expect to see 
corrections of the order of {v/cf' ~ (p/m)^ ~ 4%. If this 
is 4% of the binding energy, then it amounts to about 
300 keV. However, if relativity has a greater effect on the 
attractive a exchange part of the force (as it does in nu- 
clear matter calculations using the Walecka model) then 
we might obtain an effect 10 times larger. 

The importance of this problem has been recognized, 
and relativistic effects have been estimated using a sepa- 
rable kernel in the Bethe-Salpeter equation [Q, assuming 
minimal relativity in the Blankenbecler-Sugar equation 
and by adding corrections of first order in [v/cY to 
the Schrodinger equation All of these calculations 
include some contributions coming from relativistic kine- 
matics, but none treats the full Dirac structure of the 
nucleons, or investigates effects which might arise from a 
realistic relativistic treatment of the NN dynamics. 

The purpose of this letter is to present the first numer- 
ical calculations of the triton binding energy obtained 
from the manifestly covariant three-body spectator (or 
Gross) equations for three identical spin 1/2 particles, 
and to discuss the implications of these calculations. 
Some preliminary results were reported in conference pro- 
ceedings 0. 

The three-body spectator equations were first intro- 
duced and applied to scalar particles in 1982 and 
then extended to the case of three spin 1/2 particles in 
lectures given at the University of Hannover soon after- 
ward. Recently, a more tractable form for the equations 
has been developed, and a full derivation of the equations 
will be published elsewhere |^ . In this letter we describe 
only a few of their features briefly. 

In the absence of 3BFs the three-body scattering am- 
plitude is obtained from a sum of all successive two-body 
scatterings. Because the three particles are identical, 
each two-body scattering differs from the others only by 
a permutation, and they can therefore all be summed by 
one operator equation of the form 
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where |r-'^) is a vertex function describing the contribu- 
tion to the bound state from all processes in which the 23 
pair was the last to interact (with particle 1 a spectator), 
the two-body amplitude describes the scattering of 
the 23 pair, is the propagator for the 23 pair, and P12 
is a permutation operator interchanging particles 1 and 
2. (The factor of 2 comes from the contribution of P13 
which equals the one of P12.) 

The three-body spectator equations have the same 
structure as (|^), but incorporate the additional feature 
that the spectator is restricted to its positive energy 
mass-shell in all intermediate states. With the conven- 
tions implied above, consistency also requires that par- 
ticle 2 be on-shell, so that two particles are always on- 
shell. We think of these constraints as a reorganization 
of Eq. (^ which will, in some cases, improve its conver- 
gence. The constraints are manifestly covariant, and lead 
to the following equation 
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where the lower index labels the second on-shell particle. 
Hence only particle 3, the (unique) particle which has 
just left one interaction and is about to enter another 
one, is off-shell in Eq. (^. 

To reduce Eq. (|^) to a practical form, we take ma- 
trix elements of the operators using three-particle helic- 
ity states similar to those defined by Wick [^o). Both p- 
spin states (where p ~ + is the u spinor positive energy 
state and p = — is the v spinor negative energy state) 
of the off-shell particle must be treated. The three-body 
states will be written in the abbreviated form | Jl(23)p), 
where J is the total angular momentum of the state, p 
the p-spin of the off shell particle, 1 — {q,Xi] (where q 
and Ai are the magnitude of the three momentum and 
the helicity of the spectator in the three-body cm.), and 
(23) — {p, j, rrij, A2, A3} (where p is the magnitude of the 
relative three momentum of the 23 system, j and rrij are 
the angular momentum of the pair and its projection in 
the direction of q, and A2 and A3 are the helicities of par- 
ticles 2 and 3, all defined in the rest frame of the 23 pair). 
We will suppress all isospin indices. Using this notation, 
the final form of the three-body spectator equation for 

ri is 



(Ji(23)p|ri) = ^ ^ 

a: AiA' p' 



q dq— I dxsiTLX 
^9' Jo 



x{m)p\M'\j{2"i")p")—gP iq,p") 

tjpi' 

II I 111 

xpp^ P [1(2"3"), l'(2'3')] ^ (J'l'(2'3')p'|ri) , (3) 

where T'f/' [1(2"3"), l'(2'3')] is the matrix element of the 
permutation operator, given below, and g''{q,p) the prop- 
agator of the off-shell particle in different p-spin states 
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Because four-momentum is conserved in the relativistic 
formalism, the mass Wq of the 23 pair depends on g, 
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with Eq — \Jm^ + q2. Note that Eq. (||) includes a 
sum over intermediate helicities and angular momentum 
quantum numbers, and an integration over the internal 
spectator momentum g' and the angle x between the di- 
rections of q' and q. The integration over q' has been 
limited to the finite interval [0, ^crit], where gcrit is the 
root of the equation Wq^^.^^ — 0. At this critical spec- 
tator momentum (equal to ~ 4m/3 ~ 1200 MeV), the 
two-body subsystem is recoiling at the speed of light and 
the relativistic effects are enormous! Contributions for 

> (Zcrit are very small, and come from two-body states 
with spacelike four-momenta. It seems sensible to sim- 
ply neglect the region q' > qcHt and set the three-body 
amplitudes to zero there. As it turns out, the solutions 
go smoothly to zero as g — > gcrit anyway, so we may im- 
pose the condition that they are zero beyond this point 
without making the amplitudes discontinuous in q. 

Exchanging particles 1 and 2 implies that particle 2 
becomes the spectator and now its momentum and helic- 
ity must be expressed in the cm. frame of the three-body 
system, while the variables of particles 1 and 3 must be 
expressed in the rest frame of the 13 pair. Boosting from 
one frame to another introduces Wigner rotations of both 
the single particle and two-body helicities. The final re- 
sult for the permutation operator is 
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where the functions dmum2{l3) are the Wigner rota- 
tion matrices, and Afy,y {q,q' ,x) describes exactly the 
Wigner rotations of the off-shell particle 3, as well as 
the nontrivial matrix elements between the different p- 
spinors u and v of particle 3 as they appear in the rest 
frames of the 23 pair and the 13 pair. 

We have solved Eq. numerically for a variety of 
realistic NN models. The two-body amplitudes obtained 
for all of these models result from an exact solution of the 
two-body spectator equation, as described in Rcf. [ pT] , 
and are therefore fully consistent with the three-body 
equations. 

These models will be described in detail elsewhere. 
Briefly, they are all one-boson exchange (OBE) models 
with a kernel composed of the exchange of 6 commonly 
used bosons: the tt, ry, tr, (5, and p. The parameters of 
each model were determined by fitting to the NN phase 
shifts below 350 MeV and to deuteron properties. 



In all cases the following pion coupling was used: 
L Im 
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where p and p' are the four momenta of the incoming and 
outgoing nucleons, and the couplings proportional to 1/^^ 
do not contribute if the nucleons are on-shell. In this 
family, we fixed g'^/4:Tr = 13.34 and chose i't^ — 1, giving 
the conventional pseudovector pion coupling with large 
off-shell effects. 

A particular feature of these models, and a central 
point of this letter, is that they also include phenomeno- 
logical scalar a (with / = 0) and S {I — 1) exchanges 
with off-shell scalar-nucleon-nucleon (sNN) couplings of 
the form 



L Zm 

(to — (m — ) 



Am? 



(8) 



The vertex Ks{p' ,p) given in Eq. (|g) is the most general 
form the sNN vertex can take, but as far as we know 
the off-shell scalar couplings which depend on and Ks 
have never been studied previously. The family of models 
discussed here has = and values of v varying from 
— > 2.6, where 
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= 2.60 J/. 
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We will see that these couplings proportional to ly are 
extremely important. 

The results of our calculations are summarized in 
Fig. 1. The lower panel of Fig. 1 shows how for the fits 
to the two-body data varies with ly for this family of OBE 
models. Fits were done for values of = 0, 0.5, 1.0, 1.6, 
1.8, 1.9, 2.0, 2.2, and 2.6, and the dashed curve smoothly 
interpolates these individual cases. We emphasize that 
each of these models with different values of v are realis- 
tic in the sense that for each case OBE parameters (13 in 
all) were adjusted to give the best possible fit to the NN 
data below 350 MeV. Although the 13 parameters differ 
only slightly from case to case, the models are not quite 
equivalent. The figure shows that there is a significant 
variation in the quality of the fit; the best models lie in 
the region 1.5 < J/ < 2.0 (aU with x^ < 2.45 as shown 
in the figure). This probably rules out the model with 
ly — (for example). We conclude that the introduction 
of these i^-dependent couplings significantly improves the 
fit to the two-body data and that the implicit choice of 
1^ = made in previous work is not optimal. 

The upper panel of Fig. 1 shows the variation of the 
three-body binding energy with ly. The rapid dependence 
of the binding energy on ly is rather striking. An increase 
in ly from to 1.6 changes our prediction from —6.24 to 
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FIG. 1. Triton binding energy Et for specific values of j max 
(upper panel) and for the fits to the two-body data (lower 
panel) versus the scalar meson off-shell parameter v defined in 
the text. The curves in both panels are smooth interpolations 
through the actual calculations. The lower panel also includes 
the line = 2.45 for reference. 



—8.76 MeV, and a value in good agreement with experi- 
ment vi^ould be obtained for i> ~ 1.5, still in the range of 
i^'s vi^hich give the best fit to the two-body data. 

The panel and Table 1 (for the cases with smaller i/) 
also show how the binding energy converges as the num- 
ber of three-body partial waves, characterized by the 
highest included pair angular momentum jmax, increases. 
Because of the large increase in the predicted values as 
the number of channels increases from 28 {jmax = 1) to 
52 {jmax = 2), we were concerned about the convergence 
of the three-body calculations and studied it in detail. 
We carried the calculations all the way to jmax = 6 with 
148 channels. We find that the individual contributions 
from channels with odd j tend to cancel while those from 
channels with even j are all attractive. Thus, the steps 
from even to odd jmax are small compared to those from 
odd to even jmax- From a detailed study of the individual 
contributions we estimate that the results are converged 
to about 1 keV for = and to about 5 keV for i/ — 1.6. 

We conclude that the best of the two-body models ex- 
amined so far yield a three-body binding energy from 
about —8.5 to —9.5 MeV. In subsequent work we will dis- 
play the dependence of these results on the boosts, the 
negative p-spin states, and other relativisitc effects, and 
we will study additional two-body models. Here we will 
discuss the origin and implications of the v dependence 
which we have observed. 

To understand why the binding energy is so sensitive 
to ly, we may look at the half off-shell Born amplitude 
for scalar exchange (i.e. the amplitude with one nucleon 
off-shell). For the positive p-spin sector, we have 



TABLE I. Absolute values of the triton binding energies in 
MeV. The first row is the result when only ^5*0 and '■'Si-'^Di 
positive energy channels are included. The other rows show 
results obtained when all channels with two-body angular 
momentum j < jmax are included. The total number of 
three-body channels in each case is N. 



Jmax 


N 




coupling parameter v 








0.0 


0.5 


1.0 
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6.003 


6.345 


6.850 


7.769 
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28 


5.963 


6.318 


6.812 


7.652 


2 


52 


6.180 


6.639 


7.299 


8.441 


3 


76 


6.214 


6.695 


7.393 


8.615 
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100 


6.232 


6.726 


7.452 


8.740 


5 


124 


6.233 


6.726 


7.452 


8.736 


6 


148 


6.235 


6.731 


7.461 


8.757 
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where Vg is the usual scalar potential obtained from such 
a reduction when Vs = 0, p = {W — Ep,p) is the mo- 
mentum of the off-shell particle, W is the energy of 
the two-body system in its rest frame, and we have ig- 
nored the lower components of the Dirac spinors in car- 
rying out the reduction. The effect of the Vs depen- 
dence is to multiply the scalar potential by the factor 
Cs = [I — Vs{2Ep — W)/2m]. In applications to two- 
body scattering, the i/-dependent term is a small cor- 
rection with a sign depending on the energy, but in the 
three-body bound state it is always positive. Assuming 
an average nucleon momentum of about 200 MeV gives 
roughly a 10% variation over the range of v from to 2. 
The observed variation of about 4 MeV over this range 
is explained therefore if the average strength of the a- 
exchange potential is about 40 MeV, which is the right 
order of magnitude. This shows how the large variation 
in binding energy which we observe can be explained by 
a "small" relativistic effect. 

An OBE model with off-shell couplings has a very rich 
structure. For example, consider two successive interac- 
tions of a scalar meson with a single nucleon. The vertex 
function (|^) contains the operator m—^ which is just the 
inverse of the nucleon propagator, so that it can remove 
the nucleon propagator and contract the two interaction 
vertices to a single vertex describing the emission of two 
mesons from a single point. If the two mesons emerging 
from this point are coupled to a second nucleon they gen- 
erate a triangle or a bubble diagram. These diagrams are 
two-boson exchange terms similar to those (involving pi- 
ous) which would emerge from a nonlinear sigma model. 
Alternatively, if these two mesons couple to two different 
nucleons, they generate diagrams usually associated with 
three-body forces. It is easy to generalize this result: an 
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OBE model with off-shell couplings is equivalent to an- 
other OBE model without these couplings, but with an 
additional specific family of N -boson exchange diagrams 
and N-body forces. 

We conclude with two observations. The discovery 
that off-shell scalar couplings play an important role in 
both improving the description of two-body data and in 
predicting three-body binding energies would only have 
been possible in the context of a relativistic formalism 
closely connected to an effective field theory. The fact 
that these couplings are equivalent to a strong energy 
dependence in the context of nonrelativistic theory is 
precisely the reason they could not have been discov- 
ered there; nonrelativistic potentials are supposed to be 
energy independent. In the context of an effective field 
theory, however, these are a natural and legitimate exten- 
sion of the simplest assumption about the spin structure 
of the sNN vertex. The most general sNN vertex was 
given in Eq. (^) above, and can have only three different 
spin couplings. Once the third term depending on Ks in 
is studied, all of the possibilities will have been ex- 
hausted. In this way an effective field theory is tightly 
constrained, even if some of its interactions are strongly 
energy dependent in a nonrelativistic context. 

We believe that this way of looking at dynamics may 
very well be the most significant contribution to come 
from relativistic methods. The traditional arguments 
suggesting that relativistic effects are very small refer to 
relativistic kinematics only. As Eq. illustrates, rel- 
ativistic dynamics provides a new way to study nuclei, 
even at low energies. 
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